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We review in detail the derivation of the dielectric response function of a noninteracting system of 
spin-1/2 fermions in the random-phase approximation. Results for the response function of a Fermi 
gas in one, two and three dimensions can be obtained in closed form, and represent the baseline 
for developing a pedagogical understanding of the effect of correlations on the response functions in 
interacting systems of fermions. 



I. INTRODUCTION 

The study of quantum systems with low spatial di- 
mensions represents an active field of study that is being 
pursued aggressively in the context of emergent phenom- 
ena and has important implications for both fundamen- 
tal and applied physics. For instance, the reduced di- 
mensionality of the system is relevant to the properties 
of thin films and surfaces [ ], and for the understanding 
of high-temperature cuprate superconductors [2]. This 
in turn poses classroom challenges, aiming at exposing 
students to the array of formal techniques required to 
study such systems, with a special emphasize on linear 
response theory aimed at studying the response of the 
quantum system to a perturbation [3]. 

In quantum many-body systems the role played by 
correlations and quantum fluctuations is greatly affected 
by the reduced dimensionality of the system. Under- 
standing the role played by collective excitations in cal- 
culating the response to a perturbation of a many-body 
system represents a fundamental problem in many-body 
theory [4, 5]. The special case of a homogeneous gas 
of interacting spin-1/2 fermions has important applica- 
tions ranging from condensed matter and ultracold atom 
gases to nuclear and high-energy physics. Furthermore, 
the response of an interacting system can be represented 
as the response of a noninteracting system of particles to 
an effective self-consistent field. This is true in particular 
for weakly-interacting systems of particles. Therefore, a 
classic problem in many-body theory is the calculation 
of the response function for a system of noninteracting 
spin-1/2 fermions (e.g. electrons). 

The intent of this paper is pedagogical. The response 
function of a Fermi gas in three, two and one dimen- 
sions can be obtained in closed form, and this deriva- 
tion represents a useful classroom exercise. In addi- 
tion, the Fermi gas results are intended as the base- 
line for studying the dielectric response function in a d- 
dimensional Penn-model semiconductor [(i], perhaps the 
simplest semi-realistic demonstration of changes intro- 
duced by the presence of a gap in the single-particle en- 
ergy dispersion relation [ ]. The latter will be discussed 
elsewhere. 
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II. DIELECTRIC FUNCTION IN RPA 

Consider the interaction of an electrostatic field, 
(jjcxtj-j. i^-j^ with a system of fermions. Then, the exter- 
nal potential will induce a charge distribution, p'"'^(r,a;), 
which in turn produces a potential $'"'^(r, cj). The to- 
tal electrostatic potential $*°*(r, w) is given by the sum 
of the external potential and that associated with the 
induced charge distribution, i.e. 

$*°*(r,a;) = $°''*(r,a;) $'"^(r,tj) . (2.1) 

Similarly, the total electric charge, /9*°*(r,a;), is given by 
the superposition of the charge distribution of the non- 
interacting target, p^r,uj), and the induced charge distri- 
bution, i.e. 

p*°*(r,c.) = p(r,c.)-fp-d(r,a;). (2.2) 

The individual external, induced and total charge distri- 
butions and electrostatic potentials obey Poisson equa- 
tions of the form 

-V2$ = 47r/9. (2.3) 
In a linear-response approximation, we have 

p'"'^(r,a;) = e^J dV xir,r',io) $*°'(r',a;) , (2.4) 

where the electron response function is calculated in the 
random-phase approximation (RPA), as 

X(r,r',a;)= lim V— ^^^^^Z" ^^.S) 

r;-*0+ ^—^ri(LO — UJij) + IT] 

X V*(r)^*(r')V',(r)V'.(r'), 

or, in momentum space, 

J7-S.0+ ^ n[uj — UJij) — irj 

X (iW^-''\j){]\e-'^-'"\i) . 

Here, Si and s-^e single-particle energies and wave 

functions, and fi are occupation functions of the single- 
particle states. We have also introduced the notation. 
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hujij = Ei — Sj, and the matrix elements are calculated 



(zle'i-U) = / d'^r^Kr) e"i"-^,(r). (2.7) 



For uniform systems, we have 

$-d(q,u;) = e2v,(g) J ^ x(q,q';^) $*°*(q',^), 

(2.8) 

and Vd{q) is the Fourier transform of the coulomb poten- 
tial in d dimensions, i.e. 

r .iq.x r -27B + ln(l/(z2), if 

Vd{q) = d^x- = <^ 2TT/q , if d = 2 , 

[ 47r/q2 , if d = 3 , 

(2.9) 

with the Euler constant, 7^ — 0.577216 • • • . Therefore, 
the dielectric function in RPA reads 

e(q, q'; uj) = 5q,q- + Vd{q) x(q, c^-^)- (2-10) 

III. NONINTERACTING (FERMI GAS) MODEL 

A system of noninteracting spin- 1/2 fermions is com- 
monly referred to as a Fermi Gas (FG) system. In this 
system, one-particle states are plane waves characterized 
by momentum k, spin a ~ ±1, corresponding to the 
usual spin-up (f) and spin-down (^) fermions, and one- 
particle energies, Ck = fi.^fc^/(2m). All states below mo- 
mentum Fermi, kp, are occupied. Hence, the occupation 
of the one-particle states is given by 



nit. = Q{kF — k) 



(3.1) 



In the following the sum over states k is understood as 
an integral over all one-particle states, i.e 



d'^k 



(3.2) 



The Fermi gas density in d dimensions is calculated as 
follows 



iVi = 2 
N2 = 2 
N3^2 



k<k 



dk 


2kF 


(3.3) 






d^k 




(3.4) 


. (2^)2 


~ 2tt ' 


d^k 


7,3 


(3.5) 


. (2^)3 


" 3^ ■ 



As a matter of convenience we express momenta in terms 
of the Fermi momentum, kp, and energies in terms of the 
Fermi energy, ep = h'^kp/{2m). 



Assuming Fermi-gas-like occupation functions, /k 
nk7 the response function reads 



x(q,q;a.)-2 J ^ nk. 



(3.6) 



X [d'^k, (l-TikJ lim 



J'(k,,kj-;q,q') 

+ h{uj — ujij) + i?7 



where 9{x) is the Hcaviside function and J^(kj, k^; q) de- 
pends on the ground-state model. The factor 2 is the 
degeneracy factor corresponding to a system of spin- 1/2 
fermions. 

The matrix elements (2.7) corresponding to the case of 
a Fermi gas system are given by 



Therefore, we find that 



<5'*(q + k,:-k,) 



(3.7) 



J^(k„k,;q,q') = J^,,^ (<,k,+q " <.,k,-q) , (3-8) 

and the response (Lindhard) function for the case of the 
Fermi gas is local. We obtain 



2 f d'^k 

X0(q,^) = ^ Jj^ «k (1 - «k+q) 



(3.9) 



1 



1 



X 

Lcj - Wkq + 17] UJ + OJkq — 1»/ 

where we have introduced the notation 

Swkq = e|k+q| ^ ^ 2m^'^^ ^'^^'^ 



(3.10) 



IV. STATIC RESPONSE FUNCTIONS 

The static case corresponds to the zero-energy trans- 
fer limit, a; = 0. This case was discussed in detail by 
Kittel [ ] and we will follow his approach here. 

For w = in Eq. (3.9), it is conventional to introduce 
the notation F{q) ~ — x(q> = 0). We have: 



Fiq) = 2 

= 2 
= 2 



d'^fc 

k<k,W? 

d'^k 



k<ki 



{2TT)d 



»k(l - ni^+q) _ (1 - nk-q)nk 

Ek+q ^ Ck fik ^ Ek-q 

"k(l - "^k+q) (1 - »^k)»^k-|-q 



Ek+q — Ck 



Ck+q — Ek 



d'^fc Hk - ?lk+q 



k<k, i^Trr 



Ck+q — ek 



(4.1) 



where we changed variables in the second term based on 
the fact that the integral is over all states k. By the same 
argument, we can write 



Fici) = 2 [ ( 



/_n^ ^k 

k<kF ^^k+q — Ek Ck — Ek-q' 



(4.2) 
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whereas for q < 2, we obtain 

N2 4 kdk 



F2{q < 2) = 



N2 

(F q Jo \/q^ - 4fc2 CF 
Hence, we obtain: 

iV2 



1 - e(q-2) v/1 - (2/<z)^ 



d = 1 case: We have 



FIG. 1: (Color online) Momentum dependence of the d- 
dimensional static response functions, -F(g) = —xisi^'^ ~ 0)i 
corresponding to zero energy transfer, cj = 0. Here, momenta 
are expressed in rescaled units, i.e. in units of the Fermi mo- 
mentum, hp- 



dfc 



which gives 



2ki 

ttcf j_i q^ — 4fc2 ' 
1 + 



N^ 1 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



or 
2m 



/• d'^fc / 1 



2kq 



(4.3) 

Using rescaled coordinates, i.e. by letting k — > k/kp, we 
obtain 



In Fig. 1, we depict the momentum dependence of the 
above static response functions, Fi{q), F2{q), and ^3(5). 



V. LINDHARD FUNCTIONS 



Fici) = 



2kp 



d'^k 



2q^ 



To calculate the Lindhard function, i.e. 



ep Jk<i (27r)'^ g4_4(k.q) 



(4.4) 



We will discuss now the three-dimensional realizations of 
the static response function corresponding to d = 1, 2 
and 3. We have: 



d-^fc 



^ Jk<kp (27r) 



nk(l - Mk+q) 



(5.1) 



\UJ — UJu. 



1 



• d = 3 case: We have 
Fsiq) 



we use the formal identity 



1,3 

hip 



CF 2g 7o 



1 1-1 

dfc y 




d/i 




1 




- 4fc2^2 


UJ ±iri 


/ fcdfc log 
'0 


9- 


-2fc 
- 2fc 


•> 


Therefore, we have 



V - T i7r(5(w) 



(5.2) 



which gives 
Fz{q) 



Xlq,'^) = ^ex(q,w) + iTTOx(q,w) , (5.3) 



7V3 3 

EF 4 



log 



with 



l-ig 



• d = 2 case: We have 
»i 



F2{q) 



kdk 



d(/) 



For (7 > 2, the above gives 



(4.5) 



(4.6) 



2 f d'^k 



2u;kq 



1/ ? 2 



and 

Xmx(q,w) 



27r /■ d'^fc 

Try (2^"''('^"'^+'') 



(5.4) 
(5.5) 



F2iq > 2) = 



A^2 4 r fcdfc 
fF 9 Jo v/?^ - 4fc2 
N2 



(5(W-Wkq) +<5(CJ+Wkq) 



(1 - vi - iVir) , (4.7) 



where we introduce the notation 



^kq = Ck+q ^^k = ^('? +2k-q 



(5.6) 
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A. Real part of the Lindhard function 

In Eq. (5.4), we observe that the product nk(l — ^tk+q) 
is symmetric under the interchange k -f-)- k + q, whereas 
Wkq is odd. Hence the second term in Eq. (5.4) vanishes 
and we obtain 



^ . ^ 2 ^ /■ dU 



d-^fc 2a;kq 



(5.7) 



Using rescaled coordinates and after introducing the no- 
tation, V = hw/ep, we obtain 



72.6 x(q, 



(5., 



fc<i 



(27r)'' Vg2 _2k-q ^ _^ 2k • q 



with Qj. = v ± q^. We can also write: 

Tie x(q, J^) = Tie x'"^ (q, ~ He x'+^ (q, v) , 
where the two contributions are 

2fci „ /■ d'^/fc 1 

q| ± q • k ■ 



Jfe<i (27r) 



We note the relationship: 

7^ex(q,^^ = 0) = -^^(q) 



(5.9) 



(5.10) 



(5.11) 



We specialize now to the three-dimensional realizations 
corresponding to d 1, 2 and 3. We have: 



• c? = 3 case: We have: 



(2^)2eF Jo 



d/c 



d/i 



]^ (7^ ± 2gfc/x 



which gives 

7^ex3'^^(g, J^) 
^ TVa 3 
ep 4 



/ /cdfc log 


9| 


+ 2g/c 


/o 


9^ 


- 2qk 



(5.12) 



(5.13) 



- q| 4g2 - g,4 



2g2 



i + <z±/(2g) 



1-4/(2?) 



Hence, we obtain the well-known 3D Lindhard func- 
tion [ 1: 



'Hexaiq,^) 



TVa 3 
' ep 4 
4^2 



4g2-g,i 

-1+ o... log 



l + 9-/(2g) 



l-e/(2<7) 



log 



1 + g+/(2g) 
1-4/(29) 



(5.14) 



• d = 2 case: We have 
7^ex2'^^(g, i^) 



27r2 



l2 pi I-2TT 

^ vl kdk 



<^F Jo 



q±± 'iqk cos ( 



(5.15) 



4> 




1 — 



i 'i 

_! 5 1 



\;\7. .' 

■; 1/!/ 
;. 1/ ■ 
s I 




FIG. 2: (Color online) Momentum dependence of the real 
part of the d-dimensional Linhard functions, TZex(q,v), cor- 
responding to several values of the energy transfer, v. Here, 
momenta and energies are expressed in rescaled units, i.e. 
in units of the Fermi momentum, kp, Fermi energy, ep = 
h^kl/{2m). 



For \q%\/{2q) > 1, the above gives 



No 1 



-Re x^^\q,v) ^ '-^ - sgn(4 ) 

^ ^ ^" V^|g|/(2g)|'-fc2 



kdk 



N2 



sgn(4 



4 



2g2 



V|4/(29)|'-l 



, (5.16) 



whereas for \qW < 2q, we obtain 



TZexi^\q,>^)^ — -sgn{ql) / 

cf g Jo 



/cd/c 



4/(2g)r-fc2 



2g2 ■ 

Hence, we obtain the 2D Lindhard function 
N2 



'Hex2{q,v) = 



-1 



(5.17) 



(5.18) 



-sgn(g2)e[|4|/(29)-l] 
+ sgn(4)e[|4|/(2g)-l] 



g2 2 



2q 



- 1 



2q 



o 0.2 

ST 

-~ 0.0 



S -0.2 

^ -0.6 
-0.8 



FIG. 3: (Color online) Energy dependence of the real part of 
the d-dimensional Linhard functions, TZex{l,i'), correspond- 
ing to several values of the momentum transfer, q. Momenta 
and energies are expressed in rescaled units. 



• d = 1 case: We have: 



Ak 



■p 



Ni 1 , 



1 + 91/(29) 



i-4/(2<z) 



(5.19) 



The above gives 
log 



(5.20) 



A^i 1 



l + q^_/{2q) 



l-e/(2g) 



- log 



l + 9+/(2'?) 



l-ql/i2q) 



In Fig. 2 we depict the momentum dependence of the real 
part of the Linhard functions, 7?,ex(g, i^), corresponding 
to several values of the energy transfer, v. Similarly, in 
Fig. 3 we illustrate the energy dependence of the real 
part of the Linhard functions, TZex{q,i^), corresponding 
to several values of the momentum transfer, q. 



B. Imaginary part of the Lindhard function 

We specialize here to the case of positive energy trans- 
fer, a; > 0. Hence, the imaginary part of the Linhard 



FIG. 4: (Color online) Kinematic constraints in the energy vs 
momentum phase space for the calculation of the imaginary 
part of the d-dimensional Linhard functions. Only positive 
values of the energy, v, are physical. Momenta and energies 
are expressed in rescaled units. 



function is given by 

Imx(q,c^) (5.21) 

/ 7!r4i0(|k + q|-fcF)'^(w-Wkq). 
J<kp l^'^'j 



27r 



Using rescaled coordinates and after introducing the no- 
tation, v = hhj/ep, we obtain 



Tm x(q, V 
d'^k 



(2^)4 



(5.22) 



X |^^^e(|k + q|-l)<5(^-2k.q-g2). 

We specialize now to the three-dimensional realizations 
corresponding to d = 1, 2 and 3. We have: 



• (i = 3 case: We have: 



k'p 



k^ dk 



(5.23) 



(27r)eF Jo 

J dfie{\k + q\-l)6{i^-2kqfi-q^ 



We perform the change of variables: fi ^ p — |k + q|, 
and obtain 



Imx3{q,v) 



ep 2q 

X I kdk pdpQ{p-l)5{v-p'^ + k'^) 
J\k-q\ 



(5.24) 



Using the identity 

'5[/(p)] = ^[^(^ , with/(po) = 0, (5.25) 



6 



we obtain 



a) q = 0.5 



(5.26) 



kdk 



k+q 



dp e(p- i)S{p-po) 



\k-q\ 



where po = \/i7+W. The kinematics conditions re- 
quire that i) Po > 1, i.e. 



k > ka = \f\--v , 
and ii) |fc — g| < Po < k + q, i.e. 

-2kq < ly < + 2kq , 



(5.27) 



(5.28) 



and we recaU that v > 0. This gives the energy con- 
strain 

q'^ -2q < ly < q'^ + 2q . (5.29) 
The right-hand side of Eq. (5.28) gives the constraint 



1 /'ly 



2\q 



whereas the left-hand side of Eq. (5.28) gives 



1 fiy 

k>ki = -kr ^ ^2\q^ ^ 



(5.30) 



(5.31) 



In order to calculate Xmxaiq,'^) we need to find the 
minimum value of k, fcmin = max(fca,fcb), where kj, = 
^ I ^ — q\- Correspondingly, we introduce the notations 



^'niX3{q,v) = - 



N3 Stt 

N3 Stt ly 
e7 T 2^ 



kdk 



(5.32) 



and 



Im xliq,^) 



N3 Sir 
CF 4g jfc^ 

7V3 37r 1 I /u 



ep 4 2g 



kdk 
1 



(5.33) 



4 




FIG. 5: (Color online) Energy dependence of the imaginary 
part of the d-dimensional Linhard functions, Tmx{q,iy), cor- 
responding to several values of the momentum transfer, q. 
Momenta and energies are expressed in rescaled units. For 
convenience, we plot ^Tmxi{si,v) in the one-dimensional 
case. 



We begin by considering the case fcf, > ka- Using 
Eqs. (5.27) and (5.30), we obtain the condition 



(5.34) 



v"^ + 2q^v + g'' - 4g2 > , 

which in turn is satisfied if 

v>-q^ + 2q. (5.35) 

Hence, the relative size of ka and fcf, is obtained by com- 
paring Eqs. (5.28) and (5.35), as illustrated in Fig. 4. 



For q > 2, the inequality in Eq. (5.34) is always satis- 
fied because ly > > —q^ + 2q. Therefore, for q > 2, 



we have k„ 



kh 



1 1 E 

2 I q 



q\. We have 



Xm X3 (g > 2, 9' - 2g < < 9' + 2q) = I™ X3 (9, 



(5.36) 



For q < 2, the condition fc^ > ka is satisfied if the 
energy satisfies the condition v > —q^ + 2q. Then, we 
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obtain 



and 



Im X3{q <2,-q^ + 2q<v <q^ + 2q) = Im x§(g, v) . 

(5.37) 

If instead we have < v < —q'^ + 2q, then we have, 
ka > fcf,, and the minimum value of k is given as fcmin — 
ka — Vl ~ i^- We obtain 

Zmx3(g < 2, < 1/ < -q^ + 2q) = Imxg(g, i^) ■ (5.38) 



d = 2 case: We have: 



X / kdk 
'o "'0 



k'p 



(5.39) 



(2^)ej. 

10 e(|k + q| - I)(5(z/-2fc9cos0-g2) 



Similarly to the 3D case, we perform the change of 
variables: (f) ^ p = |k + q|, and obtain 



Imx2{q,v) = 



N2 



(5.40) 



X 4 / /cd/c 

^0 [[k + qf - p2] [p2 -[k- qf] 



Using the identity (5.25), we obtain 



Imx2{q,v) 



r kdk 

rk+q 

X / dpQ{p-l)5(p-po) 

J\k-q\ 



(5.41) 



where po — ^ v ~V fc^. The kinematics conditions and 
the p integral result in a modified range of allowed val- 
ues for /l, i.e. fcmin 

< fc < 1. The discussion regarding 
kinematics and the value of fcmin is the same as in the 
3D case. Hence, we introduce the notations 



Xmx2{q.^) = 2 



kdk 



iV2 1 ' 



(5.42) 



'1 



4 



and 



= ^ 2 



1 - 



kdk 



2 " V 



(5.43) 



Using the same arguments as in the 3D case, we find 
ItoX2('? < 2,0 < 1/ < -q^ + 2q)^Xmxl{q.v) , (5.44) 



ImX2{q>2,q'^-2q<v<q'^ + 2q)= (5.45) 
Imx2{q<2,-q^ + 2q<iy <q^ + 2q) = ImxUq, ^) ■ 



• d — 1 case: We have: 

Imxi{q,v)^~— [ dke{k + q-l)S{iy~2kq~q^) . 

EF J-l 

(5.46) 

Using the identity (5.25), we obtain 

Imxi{q,v) = ~— ^ I dk e{ko + q - 1) S{k - ko) , 
ef 4q 

(5.47) 

where ko — -^(^ — qj . The kinematics conditions re- 



quire that i) fco > 1 — q, i.e. 



>-(<?'- 2g) , 



and ii) — 1 < fco < 1, i.e. 



q -2q<iy <q^ + 2q. 



(5.48) 



(5.49) 



and we recall that i' > 0. For an energy, satisfying 
the conditions (5.48) and (5.49), we obtain 



Imxi{q,v) 



Nl TT 

cf 4g 



(5.50) 



We find that the allowed energy domains are q^ — 2q < 
V < q^ + 2q ior q> 2, and - (g^ - 2q) < v < q"^ + 2q 
for < g < 2. 

In Fig. 5 we depict the energy dependence of the 
imaginary part of the c?-dimensional Linhard functions, 
Tmx{q,i^), corresponding to several values of the mo- 
mentum transfer, q. 
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Appendix A: Useful integrals 

For the calculation of the FG response function in one 
and three dimensions we use the integrals: 



dk 



I q± ak a 



— — tanh — — — log 



1 

a 



dk 



I q^ — a^k^ aq 



= — tanh — = — log 



aq 



q + a 



q — a 



(Al) 



(A2) 
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1 g2 

/ kdk log \q±ak\= - - + log |g'| 
Jo 



4 2a2 



(A3) 



and 



kdk tanh — = — 
q 2a 



- q^ , 
1 + log 



q + a 
q — a 



(A4) 



The following integrals are used for the calculation of 
the FG response function in two dimensions: 



/.27r 

Jo 



^ ^|2^sgn(g)(M2-/c2) \ \q\>k, 

q-kcos(j) \0, \q\ < k , 

(A5) 



27T 



q"^ — k"^ cos^ ■ 



and 



kdk 



277 



0, 



^Jlf^^^^ 



1 _ ^1 _ («/<;)2 ^ (A7) 



fcdfc 



Q;2 



(A6) 



(A8) 
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